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Abstract. Let E be a flat surface of genus g with cone type singular- 
ities. Given a bipartite graph F isoradially embedded in E, we define 
discrete analogs of the 2^® Dirac operators on E. These discrete objects 
are then shown to converge to the continuous ones, in some appropriate 
sense. Finally, we obtain necessary and sufficient conditions on the pair 
F C E for these discrete Dirac operators to be Kasteleyn matrices of the 
graph F. As a consequence, if these conditions are met, the partition 
function of the dimer model on F can be explicitly written as an alter- 
nating sum of the determinants of these 2^® discrete Dirac operators. 



1. Introduction 

A dimer covering, or perfect matching, of a graph F is a cohection of 
edges with the property that each vertex is adjacent to exactly one of these 
edges. Assigning weights to the edges of T allows to define a probability 
measure on the set of dimer coverings, and the corresponding model is called 
the dimer model on T. 

Dimer models are among the most studied in statistical mechanics. One 
of their remarkable properties is that the partition function of a dimer model 
on a graph F can be written as a linear combination of 2^^ Pfafhans (deter- 
minants in case of bipartite graphs), where g is the genus of an orientable 
surface S in which F embeds. These 2^^ matrices, called Kasteleyn ma- 
trices, are skew-symmetric matrices determined by 2^^ orientations of the 
edges of F C S, called Kasteleyn orientations. P. W. Kasteleyn himself 
proved this Pfaffian formula in the planar case [12^ [T3] together with the case 
of a square lattice embedded in the torus, and stated the general fact [H] . 
A complete combinatorial proof of this statement was first obtained much 
later by Gallucio-Loebl [10] and independently by Tesler [18], who extended 
it to non-orient able surfaces. 

In [5], we studied an explicit correspondance (first suggested by Kuper- 
berg [16]) relating spin structures on S and Kasteleyn orientations on F C S. 
We also used the identification of spin structures with quadratic forms to 
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give a geometric proof of the Pfaflfian formula, together with a geometric 
interpretation of its coefficients. 

The partition function of free fermions on a closed Riemann surface S of 
genus g is also a linear combination of 2^^ determinants of Dirac operators, 
each term corresponding to a spin structure on S [Ij. Assuming that dimer 
models are discrete analogs of free fermions, one expects - in addition to 
the known relation between Kasteleyn orientations and spin structures - a 
relation between the Kasteleyn matrix for a given Kasteleyn orientation and 
the Dirac operator associated to the corresponding spin structure. 

This is well understood in the planar case with the work of Kenyon [15j. 
For any bipartite planar graph F satisfying some geometric condition known 
as isoradiality (see below), he defined a discrete version of the Dirac operator 
which turns out to be closely related to a Kasteleyn matrix of F. In particu- 
lar, its determinant is equal to the partition function of the dimer model on 
F with critical weights. In the genus one case, the following observation was 
made by Ferdinand [9j as early as 1967: For the M x N square lattice on the 
torus with horizontal weight x and vertical weight y, the determinants of 
the four Kasteleyn matrices behave asymptotically, in the MN — )• oo limit 
with fixed ratio M/N, as a common bulk term times the four Jacobi theta 
functions 6k{0\T), where r = i^^- This reproduces exactly the dependance 
of the determinant of the Dirac operators on the different spin structures 
observed by Alvarez-Gaume, Moore and Vafa [1]. 

The higher genus case remains somewhat mysterious. The only results 
available are numerical evidences, for one specific example of a square lattice 
embedded in a genus two surface, that the determinants of the 16 Kasteleyn 
matrices have a dependance that can be expressed in terms of genus two 
theta functions [6j. 

In short, it is fair to say that relatively little is understood of the expected 
relation between Kasteleyn matrices and Dirac operators on surfaces of genus 
g > 2. With this paper, we aim at filling this gap. Here is a summary of 
our results. 

We start in Section [2] by defining a discrete analog of the d operator on 
functions on a Riemann surface S. Because there is no "canonical" such 
discretization, some geometric conditions may be naturally imposed on the 
pair F C S. Following Duffin [7], Mercat [T7], Kenyon [15J and many others, 
and for reasons that will become apparent along the way, we work with 
bipartite isoradial graphs. More precisely, we encode the complex structure 
on S by a flat metric with cone type singularities supported at S C S, and 
consider locally finite graphs F C S with bipartite structure V^(F) = BUW 
satisfying the following conditions: 
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(i) Each edge of F is a straight hne (with respect to the flat metric on 
S), and for some positive 5, each face / of F C S contains an element 
Xf at distance S of every vertex of /. 

(ii) A singularity of E is either a black vertex of F, or a vertex Xf the 
dual graph F*. 

For such a pair F C S, we introduce a discrete B operator defined on C^, 
and call / G discrete holomorphic if Bf = 0. (See Definition 12. li ) 
This operator satisfies some natural properties (Proposition 12. 2p and ex- 
tends previous constructions of Mercat [17], Kenyon p!5] and Dynnikov- 
Novikov |8j. Note that these authors impose strong conditions on F in 
order for d to be defined: F needs to be the double of a graph in [17J (and 
therefore, does not admit any perfect matching in higher genus cases, see 
Remark I2.3p . it needs to be planar in [15j . while only the triangular (or 
dually, the hexagonal) lattice is considered in [8]. On the other hand, our 
discrete B operator imposes essentially no combinatorial restriction on F: 
any locally finite bipartite graph such that each white vertex has degree at 
least three can be isoradially embedded in an orientable flat surface S with 
conical singularities S C V{T*) U B (Proposition I2.4( ). This section is con- 
cluded with a convergence theorem: if a sequence of discrete holomorphic 
functions converges to a function /: S — )• C in the appropriate sense, then 
/ is holomorphic. (See Theorem 12.51 for the precise statement.) 

In Section [3l we twist the discrete B operator on F C S by discrete spin 
structures A to obtain 2^^ discrete Dirac operators 

provided each cone angle of S is a multiple of 27r (Definition 13. 9p . The 
convergence theorem then takes the following form: let be a sequence 
of discrete spin structures on S discretizing a fixed spin structure L. If a 
sequence tpn of discrete holomorphic spinors (that is: D\^ipn = 0) converges 
to a section ■0 of the line bundle L — )• S, then -0 is a holomorphic spinor. 
(See Theorem 13.121 ) 

Section [3] contains the core of this paper. First, we extend the Kasteleyn 
formalism from {ibl}-valued flat cochains on F C S (that is, Kasteleyn 
orientations) to G- valued ones for any subgroup G of C*. We believe that 
the resulting existence statement (Proposition 14. 2( ) and generalized Pfaffian 
formula (Theorem 14.31 and Corollarv 14. 4p are of independant interest. We 
use them to prove our main result: 

Theorem. Let Ti be a compact oriented flat surface of genus g with conical 
singularities supported at S and cone angles multiples of 2-k. Fix a graph 
F with bipartite structure V{r) = B L\ W , isoradially embedded in S so 
that S C B Li V(T*). For an edge e ofT, let v{e) denote the length of the 
dual edge. Finally, let D\ : — > C'^ denote the discrete Dirac operator 
associated to the discrete spin structure A. 
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There exist 2^^ non- equivalent discrete spin structures such that the corre- 
sponding discrete Dirac operators {D\\\ give non- equivalent Kasteleyn 
matrices of the weighted graph (r, i/), if and only if the following conditions 
hold: 

(i) each conical singularity in V{T*) has angle an odd multiple of2iT; 
(a) for some (or equivalently, for any) choice of oriented simple closed 
curves {Cj} in T representing a basis of Hi{T,;Z), 

is a multiple of 27: for all j , where a^,(C) denotes the angle made by 
the oriented curve C at the vertex v as illustrated below. 




V 



As a consequence, given any graph F C S satisfying the conditions above, 
the partition function for the dimer model on (F, v) is given by 

^(r,^^) = ^| E (-1)^'''^'^ det(Z),) , 
Aes(s) 

where Arf(A) € Z2 denotes the Arf invariant of the spin structure A (Theo- 
rems 14.91 and 14. lip . Our final result - Theorem 14.141 - states that the Dirac 
operators on any closed Riemann surface can be approximated by Kasteleyn 
matrices. More precisely, for any closed Riemann surface of positive genus, 
there exist a flat surface S with cone type singularities inducing this com- 
plex structure, and an isoradially embedded bipartite graph T C S, with 
arbitrarily small radius, satisfying all the hypothesis and conditions of the 
theorem displayed above. 

Acknowledgements. The author would like to thank Mathieu Baillif, Gio- 
vanni Felder and Nicolai Reshetikhin for useful discussions. 

2. The discrete B operator on Riemann surfaces 

The aim of this section is to introduce a discrete analog of the d operator 
on functions on a Riemann surface, extending works of Duffin [7], Mer- 
cat |17j and Kenyon |15] . As this definition requires a substantial amount of 
notation and terminology, we shall proceed leisurely, starting by recalling in 
Paragraph 12 . 1 1 the main properties of flat surfaces with conical singularities. 
We then give in Paragraph 12.21 discrete analogs of all the geometric objects 
involved in the definition of d (see Table [l]) . This will lead up in Paragraph 
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12.31 to the - by then, quite natural - definition of the discrete operator. The 
section is concluded with a convergence theorem (Theorem l2.5l in Paragraph 
12. 4p . justifying further our definition. 

2.1. Flat surfaces. Our discrete d operator will be defined for graphs em- 
bedded in so-called flat surfaces with conical singularities. Since these ob- 
jects are ubiquitous in the present paper, we devote this first paragraph to 
their main properties, referring to [19j for further details. 

Given a positive real number 0, the space 

Ce = {{r,t) : r > 0, t G M/0Z}/(O,t) ~ (0,t') 

endowed with the metric ds^ = dr^ + r'^dt^ is called the standard cone 
of angle 6. Note that the cone without its tip is locally isometric to the 
Euclidean plane. Let S be a surface with a discrete subset S. A flat metric 
on S with conical singularities of angles {Ox}x€S supported at S is an 
atlas {(j)x : Ux ^ U'^ C Ce^^x&s-, where Ux is an open neighborhood of x G 5, 
(j)x maps X to the tip of the cone Cq^ , and the transition maps are Euclidean 
isometries. 

This seemingly technical definition should not hide the fact that these 
objects are extremely natural: any such flat surface can be obtained by 
gluing polygons embedded in along pairs of sides of equal length. For 
example, a rectangle with opposite sides identified will define a flat torus 
with no singularity. On the other hand, a regular Ag-gon with opposite sides 
identified gives a flat surface of genus g with a single singularity of angle 
2n{2g — 1). In general, the topology of the surface is related to the cone 
angles by the following Gauss-Bonnet Formula: if S is a closed flat surface 
with cone angles {Ox}xeSi then 

J^(27r-e,.) =2^X(S). 

xes 

For the purpose of this paper, the most important property of fiat metrics 
is that they encode complex structures on oriented surfaces. Indeed, the 
conformal structure on S \ 5" given by a flat metric extends to the whole 
oriented surface, deflning a complex structure on S. Furthermore, let S be a 
closed oriented surface, 5 C S a discrete subset, and {6x}xes a set of positive 
numbers satisfying the Gauss-Bonnet Formula. Then, for each complex 
structure on S, there exists a flat metric on S with conical singularities of 
angles {9x}x£S supported at S inducing this complex structure. 

For example, any complex structure on the torus can be realized by a 
flat surface with no singularity: simply consider the parallelogram in the 
complex plane spanned by the pair of periods of the torus, and identify 
the opposite sides. Similarly, deforming the regular octagon in such a way 
that the sides are organized into pairs of equal length allows to realize any 
complex structure on the genus two surface. 
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2.2. Some discrete geometry. Let us begin by briefly recaUing the defini- 
tion of the d and B operators on a Riemann surface S. (This will also fix some 
notation). The complex structure J on S induces a decomposition of the 
complexified tangent bundle TS*^ into TS"*" ©TS~, and therefore, a decom- 
position of complex- valued vector fields C^{TTF) = C°°(TS+)eC°°(TS-). 
The elements of C°°(rS+) (resp. C°°(TS")) are the vector fields for which 
the action of J is given by multiplication by i (resp. — i). Similarly, the com- 
plex cotangent bundle splits, resulting in a decomposition of the complex- 
valued 1-forms on E: 

Q^(s, c) = c°°{t*yF) = c°°(r*E+) e c°°(t*e-) = e j^°'^(e). 

Note that the forms of type (1,0) (resp. (0,1)) are the 1-forms such 
that for any vector field V, ip{J{V)) is equal to i^iV) (resp. —i(p{y)). 
Finally, the exterior derivative d: C°°(E) — ?• r2^(E) induces a C- linear map 
d^: C'^(E,C) ^ ^]i(E,C) = f^i'0(E)®JlO'i(E), whose composition with the 
natural projections defines the Dolbeault operators d: C°°(E,C) r2^''^(E) 
and d: C°°(S,C) ^ OOji(S). Recall that a function / G C°°(S,C) is holo- 
morphic if and only if df is zero. 

We are now ready to start our discretization procedure. First and fore- 
most, a Riemann surface E is a surface. To encode the topology of E, fix 
a locally finite graph P C E with vertex set ^(r) and edge set -E(r), such 
that E \ F consists of disjoint open discs. In other words, F is the 1-skeleton 
of a cellular decomposition of E. For notational simplicity, we shall assume 
throughout this section that F has neither multiple edges, nor valency one 
vertices. (Note however that all our results hold in the general case as well.) 

As explained in the previous paragraph, a standard and beautiful way to 
encode a complex structure on an oriented surface S is to endow this surface 
with a flat metric with conical singularities. Note that any point in E \ S' 
has a well-defined tangent space. In particular, the space X(E) of vector 
fields on E can be naively discretized by X(D), the space of vector fields 
along some discrete subset C E \ 5. 

Next, we wish to encode in the graph F the decomposition of j£(E)'^ = 
C°°(rE'^) induced by the almost complex structure. A convenient way to 
do so is to consider a decomposition ViV) = B UW oi the vertices of F 
into, say, black and white vertices, together with a perfect matching M on 
F pairing each black vertex b E B with a white one w eW and vice versa. 
Any perfect matching will do the job, provided F is bipartite, that is: no 
edge of F links two vertices of the same color. Hence, when the surface 
E is endowed with an almost complex structure, it is natural to consider 
a bipartite graph F C E together with a perfect matching M on it. The 
discrete analog of the decomposition C°°{TT,^) = C°°(rE+) © C°°{TT,-) is 
then given by X{DM f = X{B)(BX{W), where Dm C E denotes the discrete 
set consisting of the middle point of each edge in M. The complex structure 



Discrete Dirac operators and Kasteleyn matrices 

Table 1. Discretization dictionary, part 1 



7 



the geometric object 


the discrete analog 


a surface T, 


a graph F C S inducing a cellular 

QtJCUllipublLlUli (Ji Zj 


a conformal structure on S 


a flat metric on S with conical 

blllgUiai ILlCb O Zj 


the space of vector 
iicicib on Zj 


the space X{D) of vector fields along 
some Qibciexe suubex u \^ 


the decomposition 
O [11^ ) — 

induced by an almost 
complex structure on S 


a bipartite structure T^(r) = BUW 
together with a perfect matching IvI on 
the graph F, inducing a decomposition 


the space of 
(1, 0)-forms on S 




the space O'^'^(S) of 
(0, l)-forms on T, 




the space C°°(S,C) of 
complex functions on S 


~ C^, identified via M 


jJpdFdxdy = -^jQpFdz 


the definition of the discrete d operator 



on X{B) (resp. X{W)) is such that multiplication by i corresponds to the 
90 degrees rotation of the tangent vectors in the positive (resp. negative) 
direction, which we will represent in our figures as counterclockwise (resp. 
clockwise). 

In the same way, r2^''^(S) will be encoded by the space Q^{B) := HfoGB X{b)* 
and JlO'i(S) by ^^^(l^) := UweW^i^)*^ where X{v)* denotes the dual to 
the 1-dimensional complex vector space X{v) of tangent vectors at the ver- 
tex V € V(^)- Finally, the space C°°(S,C) can be discretized both by 
and by C^, which are identified via the perfect matching M. 

Table [T] summarizes our notations and the dictionary between the geo- 
metric and discrete objects considered in this paragraph. (The last entry 
will be explained shortly.) 

2.3. The discrete d operator. Let F be a bipartite graph embedded in a 
flat surface S with conical singularities supported at S. According to the 
discussion above, the operator B: C°°($],C) — )• r2'^'^(S) should discretize to 



8 



David Cimasoni 



a C-linear map 

B: ^n\w) = Yl X{w)*. 

These spaces make sense as soon as no white vertex is a singularity, so let 
us only assume S C T,\W for now. 

Following the terminology of Kenyon (1^ . we shall say that T is isora- 
dially embedded in S if each edge of F is a straight line, and if for some 
5 > 0, each face / of F C S contains an element Xf such that d{xf,v) = 5 
for all vertices v oi df. We shall furthermore assume that a singularity of S 
is either a black vertex b of F, or a vertex of the dual graph F*, that is: 
S C V{T*)UB. 




Figure 1. A isoradial graph F (black vertices, solid edges), 
its dual graph F* (lighter vertices and edges), and the asso- 
ciated rhombic lattice (all vertices, dashed edges). 

Given an isoradially embedded graph F C S, the associated rhombic 
lattice is the graph i?r with vertex set V{T)L)V{T*) and edges joining each 
vertex of F with the center of the adjacent faces, as illustrated in Figured) 
Since F induces a cellular decomposition of S, so does the rhombic lattice 
Rr C S. Furthermore, as the singularities of S lie among the vertices of 
the rhombic lattice, one easily checks that the faces of this lattice are actual 
planar rhombi. Therefore, the metric space S should be understood as 
planar (paper) rhombi pasted together along their boundary edges. 

For a fixed white vertex w G W, let St{w) C S denote the star of w in the 
rhombic lattice, that is, the union of all the closed rhombi adjacent to w (see 
Figure [2]). As w does not belong to the singular set either, the whole star 
can be isometrically embedded in the Euclidean plane, as "demonstrated" 
by cutting and pasting paper rhombi. Hence, one should really think of S 
as planar stars as in Figure [2] pasted along their edges. (Note that since 
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multiple edges and valency one vertices are not allowed, we avoid the case 
where boundary edges of a star are glued together. But again, the difficulty 
of the general case would only be notational.) 




bjn 

Figure 2. The star St{w) of the white vertex w. 

As mentioned in the introduction, isoradial graphs first appeared in the 
work of Duffin [7j (in the form of planar rhombic lattices) as a large class of 
graphs for which the Cauchy-Riemann operator admits a nice discretization. 
It should therefore not come as a surprise that such a condition is necessary 
for our definition. 

Definition 2.1. Let F be a bipartite graph isoradially embedded in a flat 
surface S with conical singularities S C V{T*) U B. Given an edge (w, b) of 
F, let I'iw, b) denote the length of the dual edge. The discrete d operator 
is the linear map 9: -> Q.^{W) = YiweW ^(^)* defined by 

(a/)H(v») =: (AJ)(v.} = ^Jlb) 

for / € C^, w ^ W and Vw G X{w). The sum is over all vertices b adjacent 
to w, and '&vi'w, b) denotes the angle at w € W from the tangent vector 
to the oriented edge {w,b), as illustrated in Figure [3l 

A function / € is discrete holomorphic if Bf = 0. 




Figure 3. Definition of the angle i3v{w,b). 
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This definition should be understood as a discretization of the formula 
// dF{x + iy) dx dy = / F{z)dz. 

J Jp 2 Jgp 

Indeed, given / G C^, let /: S ^ C be defined (almost everywhere) by 
f{p) = f[h) if p belongs to the interior of the star St{b). Fix a white vertex 
w, and consider an isometric embedding (j) of the corresponding star St{w) 
into C. Setting = iTw(l))~^ {!) ■, F = f o cp'^, and using the notations of 
Figure [21 we get 



{d^f){V^) « {dFmw)) 
1 



Area{St{w)) 
—i 



dF{x + iy) dx dy 



F(z)dz 



2Area{St{w)) J^(9st(w)) 

^ m 

j=i 

2Area{St(w)) 

This definition extends previous work of Duffin [7] , Mercat [T7] , Kenyon [15] , 
Dynnikov-Novikov [8j and Chelkak-Smirnov [3], as described below. 

Special case 1. If the fiat surface S has trivial holonomy, there is a well- 
defined constant vector field in X{W). Evaluating the discrete d operator 
at this vector field yields a map K: — >■ C^. In the special case where S 
has no singularity (which is only possible if S is the plane, a cylinder or a 
torus), this map coincides with the discrete d operator on planar bipartite 
isoradial graphs defined by Kenyon [l_5j (up to the normalization constant). 
For the planar hexagonal lattice, this operator is conjugate to the discrete 
d operator considered by Dynnikov-Novikov in [8]. 

Special case 2. Let 9 be a (non-necessarily bipartite) graph embedded in 
a flat surface S with singularities S C V{9)UV{5*) =: A. The double of 9 
is the bipartite graph r = 9u9*c5] with black vertices B = A and white 
vertices W = £'(9) n£'(9*) =: 0, or the other way around. If 9 is isoradially 
embedded in S, then so is 9* (with same radius 6) and F (with radius 5/2). 

In the special case of planar double graphs, an element / G is either 
a function on the vertices of the associated rhombic lattice (when B = A), 
or a function on the set of rhombi (when B = (}). As E = C, our discrete d 
operator can be evaluated at the constant vector field 1 € TpC = C, yielding 
two maps — > and — )■ C^. These correspond exactly to the two 
discrete d operators defined by Chelkak and Smirnov in [3]. 
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Special case 3. Finally, let us con- CR.{ps,eps} not found (or no BBox) 
sider the more general case of double 

graphs isoradially embedded in flat surfaces with conical singularities, with 
bipartite structure B = A. Let w (z W = he a fixed rhombus, and 
let x,y,x',y' denote its vertices enumerated counterclockwise, as illustrated 
above. Then, the equality d^f = coincides with the very intuitive "dis- 
crete Cauchy-Riemann equation" studied by Duffin [7] (in the planar case) 
and Mercat [17\ : 

f{y') - f{y) _ , fi^') - fi^) 

d{y,y') d{x,x') 

Let us mention several natural properties of the discrete d operator. 
Proposition 2.2. The discrete d operator satisfies the following properties. 
{'>■) If f ^ is constant, then df = 0. 

(a) Given a fixed white vertex w, let f G ([^BnSt(w) restriction of 

a coordinate chart on a neighborhood of St{w). Then, d^f = . 
(Hi) Let f G ([^Bnst^w) ^j^j^j above. Then its complex conjugate f 

satisfies {dwf){Vw) = \Vw\. 

Proof. Let / € be a constant function. Fix a white vertex w ^ W 
and a unit tangent vector € X{w). Let cp: St{w) — )• C be an isometric 
embedding mapping w to the origin and Vw to the direction of the unit 
vector 1 € C = TqC With the notation of Figure O observe that for all j, 

uiw,b,)e'''y'^^''^^ = i{cP{xj) - Hxj^i)), 

as both these complex numbers have same modulus and argument. It follows 
that 

2Area[St{w)) ^ 

proving the first claim. To check the second one, let (j): U ^ C be a co- 
ordinate chart with St{w) C U, and let Vw G X{w) he the unit tangent 
vector corresponding to the edge {w,Xm) (recall Figure [2]). By the first 
point above, it may be assumed that 4'{w) = 0. Clearly, one can also as- 
sume that Tyj(t>{yw) = 1 G TqC For j = 1, . . . m, let aj denote the angle at 
w of the rhombus corresponding to the edge {w,bj). Note that 

sin(aj) = 2sin(aj/2)cos(aj/2) = ^ ^'^ ^' . 
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Since f{bj) = (j){bj) = d{w,bj)e^^y^'"'^i\ we get 

m 

2Area{St{w)){d^f){V.^) = Y.v{w,h,)e'''^^^^''^^ f{bj) 

m 

i=i 

m 

i=i 

m 

i=i 

using the fact that YlT=i = 2vr. Finally, 

\V \ 

showing the third claim. □ 

Remark 2.3. As pointed out in the special cases above, most authors have 
considered discrete d operators defined on double graphs only. It is however 
crucial for us to consider more general graphs, for the following reason. In 
Section SI we shall turn to the problem of counting perfect matchings on 
a (finite) bipartite graph T embedded in a (compact) surface S. For such 
a matching to exist, one necessary condition is that the number of black 
vertices equals the number of white ones. But in the case of a double graph 
r = D(S) C S, this condition gives 

Q = \B\- \w\ = |y(g)| + |F(g)| - \e{^)\ = x(s). 

Hence, no double graph as above admits a perfect matching unless S is a 
torus. 

On the other hand, our setting imposes almost no restriction on the com- 
binatorial type of the graphs considered: 

Proposition 2.4. Let T be a locally finite bipartite graph such that each 
white vertex has degree at least three. Then, T can be isoradially embedded 
in an orientable flat surface S with conical singularities S C V{T*) U B. 

Proof. For each w € W, fix a cyclic ordering of the m adjacent edges (so 
that multiple edges are consecutive) and form the symmetric star St{w) by 
pasting together m rhombi of side length 6 and of angle In jm according to 
this ordering. Note that the surface St{w) is endowed with an orientation 
given by the cyclic ordering. In case of multiple edges or black vertices of 
degree 1, identify the corresponding boundary edges of St{w) accordingly. 
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(This respects the orientation of the star.) For each 6 G i?, fix a cychc 
ordering of the adjacent edges, and glue the stars St{w) along their boundary 
edges according to these orderings, in the unique way compatible with the 
orientations of the stars. The result is an oriented flat surface S with conical 
singularities supported at S. By construction, F is isoradially embedded in 
S and S is contained in V(T*) Li B. □ 



As an example, consider the complete bipartite graph 
-f^3,3- With a natural choice of the cyclic orderings 
around the vertices, the construction above yields the 
honeycomb lattice embedded in the flat torus illustrated 
opposite. (The pairs of opposite sides of the big hexagon ^ 
are identified.) 

To conclude this paragraph, note that if S C S \ i?, one can define the 
discrete d operator as the C-linear map d: — )• Q^{B) = Yli^^B^i^)* 
given by 

for g G C^, b (z B and C/;, G X(6). Here again, the sum is over all vertices 
w adjacent to 6, and ^u{w-, h) denotes the angle at 6 G -B from the tangent 
vector Uh to the oriented edge {w,b). This construction generalises the one 
given in [15], which corresponds to the case with no singularity. However, 
the discrete d operator being sufficient for our purposes, we shall not study 
d in the present paper. 

2.4. A convergence theorem. The aim of this paragraph is to prove the 
following result. 

Theorem 2.5. Let T, be a flat surface with conical singularities supported 
at S. Consider a sequence F„ of bipartite graphs isoradially embedded in S 
with S C 1^(F*) U Bn- Assume that the radii 5n ofVn converge to 0, and 
that there is some rj > such that all rhombi angles of all these F„ 's belong 
to [rj, IT — r]]. Let fn G C^" be a sequence of discrete holomorphic functions 
converging to a function f : S — )• C in the following sense: for any sequence 
Xn G Bn converging in S, the sequence fn{xn) converges to /(lim„x„) in C. 
Then, the function f is holomorphic in S. 

Our proof will follow the same lines as the one of Mercat |17l pp. 192-195], 
a notable exception being the discrete Morera Theorem below (Lemma l2.8p . 
Let us start with a straightforward generalization of [ITj Lemma 2]. 

Lemma 2.6. Let X be a metric space. Consider a sequence of functions 
fn'- X ^ C converging to f : X ^ C in the following sense: for any conver- 
gent sequence Xn in X, the sequence fn{xn) converges to /(lim„x„) in C. 
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Then, the function f is continuous, and is the uniform limit of fn on any 
compact. 

Proof. To show that / is continuous at an arbitrary point x € X, pick a 
sequence Xj converging to x in X. For any j = 1,2,..., the hypothesis 
apphed to the constant sequence Xj yields the existence of an index rij such 
that \fnj{xj) — f{xj)\ < Let y„ be the sequence given by y„ = Xj if 
n = Uj, and yn = x else. As yn converges to x, fniVn) converges to f{x) 
and so does the subsequence fnj{xj). It follows that 

- f{x)\ < \f{x,) - + |/n,(x,) - fix)\ 

is arbitarily small, proving the first claim. 

To show the second one, let us assume by contradiction that /„ does not 
converge uniformly on some fixed compact C C X. This would imply the 
existence of a convergent sequence x„ in C with \fn{xn) ~ f{xn)\ greater 
than some e > for all n. On the other hand, the hypothesis together with 
the continuity of / at x = lim„ Xn imply 

\fn{Xn) - f{Xn)\ < |/„(Xn) - f{x)\ + \f{x) - f{Xn)\ < E 

for n big enough, a contradiction. □ 

Lemma 2.7. Let V he a graph isoradially embedded in the Euclidean plane, 
such that all rhombi angles belong to the interval [r],Tr — ri] for some 77 > 0. 
Then, for any vertex v of T, and for any two elements x, x' in the boundary 
dSt{v) of the star of v, 

ddst{v){x,x') ^ 2n 
\x — x'\ ~ r]sm{r]/2) 

Proof. To simplify the notation, let S stand for the star St{v) throughout 
this proof. If X and x' belong to the same rhombus of S, then the quotient 
above is easily seen to be maximal when x and x' are at the same distance of 
the vertex opposite to v. In such a case, this quotient is equal to 1/ sin(Q/2), 
where a denotes the angle of this rhombus at v. Since 77 < a < tt — 77, it 
follows 

dds{x,x') ^ 1 ^ 1 ^ 27r 
|x — x'l ~ sin(Q;/2) ~ sin(r7/2) ~ ?7sin(77/2)' 

Let us now assume that x and x' belong to adjacent edges, but distinct 
rhombi of S. If the corresponding angles at w are equal to a and a', then 
the argument above gives the inequality 

dds{x,x') ^ 1 ^ 1 ^ 27r 

|x — x'l ~ sin((a + a')/2) ~ sin(77) ~ r/sin(?7/2) 

If X and x' lie on adjacent rhombi, but non-adjacent edges of the star, then 
|x — x'l is bounded below by Ssm{rf), where 6 denotes the length of the 
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rhombus edges. On the other hand, dQs{x,x') < 4:5 a,s x and x' belong to 
adjacent rhombi. Therefore, in this case 

dds{x,x') ^ 4 ^ 27r 



sin(r/) r/sin(r//2) 

Finally, consider the case where x and x' do not belong to adjacent rhombi. 
Fix a rhombus between them, and let a denote its angle at v. This time, 
Ix — x'l is bounded below by 2(5sin(a/2) > 25sin(r//2), while the distance in 
dS is bounded above by half of the length of 95, that is 



ddsix, x) < '- = #{rhombi m S"} • 5 < . 

2 1] 



This implies 



das{x,x') ^ TT ^ 2tt 



\x — x'\ rism{ri/2) rj s'm{ri / 2) '' 
and concludes the proof. □ 

The last lemma requires some preliminaries. As above, let F be a bipartite 
graph isoradially embedded in a flat surface S. Given a function / : S — )• C 
and a white vertex w of F, set 



/:= / /(7(t))(0o7y(t)dt, 
ast{w) Ja 

where 7: [a, b] St^w) is a parametrization of dSt{w) and (p: St{w) C is 
an isometric embedding. Obviously, the value of this integral depends on the 
choice of the chart </>. However, the choice of another chart would multiply 
the result by a modulus 1 complex number. In particular, the vanishing of 
this integral does not depend on such a choice, and the following statement 
makes sense. 

Lemma 2.8 (discrete Morera Theorem). Let T be a bipartite graph isoradi- 
ally embedded in a flat surface S with conical singularities S C V{T*) U B. 
Given f E C'^, let f : T, ^ C be the function defined by f{p) = ^ Si=i fi^j) 
if p belongs to CfjLi ^^(^j)- '^hen, f is discrete holomorphic if and only if 

Proof. Let w be a white vertex, and let (p: St{w) — > C be an isometric 
embedding of the corresponding star. Fixing a unit vector € 3i{w), one 
can assume that T^cp maps 14) to 1 E C. With the notation of Figure [H we 
get the equality 

„ m 

/ / = y,(.H^j) - <P{xj-i))f{bj) = 2iArea{St{w))id^f)iV^), 
Jast(w) jr{ 

and the lemma follows. □ 
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Proof of Theorem \2.5[ As in Lemma 12.81 extend € C^" to a function 
/„: S ^ C by setting Up) = ^ ET=i U^j) if P belongs to 07=1 St{bj). 
By assumption, there is a function / : S — )• C such that, for any sequence 
Xn G Bn converging in S, the sequence fn{xn) converges to /(hm„a:„) in C. 
We claim that this statement remains true for the extensions fn- S — )• C. 
Indeed, let us fix a sequence Xn in S converging to x. For all n, there 
exist black vertices bn\ . . . ,hrT'^ such that x„ belongs to the star St{hn^) 
for j = 1, . . . , m. Let 6„ (resp. b'^) be one of these vertices where Re fn is 
maximal (resp. minimal) on this set. By definition, 

Refn{h'^) < RefniXn) < Refnibn). 

Since 6„, 6^ and Xn all belong to the adjacent (or identical) stars St{bn) and 
Stib'n) whose diameter is at most 45ji, and since 5n converges to zero, both 
sequences 6„ and b'^ converge to x = lim„ Xn- By the assumption, fnipn) and 
fn{b'n) both converge to /(x). By the inequalities displayed above. Re fn{xn) 
converges to Re f{x). The same argument shows that Im fn{xn) converges 
to Imf{x), proving the claim. 

Lemma 12.61 asserts that / : S C is continuous and the uniform limit 
of /„: S ^ C on any compact. Since the singular set S* C S is discrete 
and / is continuous, it is now sufficient to check that / is holomorphic on 
S \ (S*. Hence, we can restrict ourselves to (simply-connected) domains of a 
Euclidean atlas for the flat surface 5 C S. In other words, we can assume 
that we are working in a simply connected planar domain U C C Finally, 
by Morera's theorem, it is enough to show that f{z) dz vanishes for any 
piecewise smooth loop 7 in U . 

So, let 7 be such a loop, let I denote its length, and let n be a fixed index. 
Each time 7 meets some star St{wn), entering it at a point x and leaving 
it at x', replace 7 n St{wn) by the path in dSt{wn) realizing the minimal 
distance in dSt{wn) between x and x' . This yields a new loop 7^ contained 
in StY^-, the union of all rhombus edges adjacent to black vertices of r„. By 
Lemma 12.71 its length satisfies 

where M(ry) stands for the uniform bound ^ sin{ri/2) • diameter of a star 

St(wn) is at most 45^, the union of all these stars meeting 7 is contained 
in the tubular neihborhood of 7 of diameter 85n, which also contains the 
compact set C„ enclosed by 7 and 7„. Therefore, Area{Cn) < 8dni- 

We shall now prove that the sequence f{z) dz converges to f(z) dz. 
Let us first assume that / is of class C^. In such a case, df is bounded 
above by some constant M on the compact set C given by the tubular 
neighborhood of 7 of diameter 8 niscKfi 6fi. As C contSiins 3>11 th.6 C^^'s, this 
yields a uniform bound for df on all CnS. By Stokes formula and the 
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inequality above, 



/ fiz)dz- [ f{z)dz 




/ f{z) dz 




1 1 df{z) dz A dz 






JdCn 







< 



\df{z)\dzAdz < 8 6niM, 



Cn 



proving the claim in this special case. In the general case of a continuous 
function /, let gu be a sequence of functions on U converging uniformly 
to / on every compact. By the inequalities displayed above, we obtain that 



I'-l ' 


< 


f - 9k 


+ 


9k- 9k 


+ 


/ 9k- j f 


J7 J7„ 




J 'J J y 




Jin 




J^n J^n 



< £ sup[/ + 
c 

is arbitrarily small, proving the claim. 



\6n£M + M{rj)isnp\f-gk\ 
c 



Recall that /„ converges uniformly to / on the compact C. Therefore, for 
any fixed index k, 



fn{z) dz 



Ik 



f{z)dz 



Ik 



<M{7])£sup\fn-f\ 

c 



is arbitarily small, so /^^ /„ converges to /^^ /. 

We are finally ready to show that f[z) dz is equal to zero. As Str^ 
induces a cellular decomposition of the simply-connected domain [/, the 
cycle 7„ C 5tr„ is a cellular boundary, that is, 7„ = 9 (^^^ St{wn)) for 
some vertices Wn- Since /„ is discrete holomorphic, Lemma 12.81 implies 

/ Uz) dz = J2 [ ?n{z) dz = 0. 
By the three claims above. 



7fc 



k n 



f{z) dz = lim / f{z) dz = limlim / fn{z) dz = lim / fn{z) dz = 0. 



7fe 



This concludes the proof of the theorem. 



□ 



3. Discrete Dirac operators on Riemann surfaces 

In the previous section, we defined a discrete analog of the d operator 
on functions on a Riemann surface S. The aim of the present section is 
to modify this construction, yielding an analog of the Dirac operator D 
on spinors on S. Here again, we shall start by giving in Paragraph 13.11 
discretizations of all the geometric objects involved in the definition of D 
(Table [2]). The actual definition of the discrete Dirac operator is to be 
found in Paragraph 13.21 while Paragraph 13.31 deals with the application of 
our convergence theorem to spinors (Theorem I3.12p . 
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3.1. More discrete geometry. Let us first recall the definition of the 
Dirac operator on a closed Riemann surface S, referring to [2j for details. Let 
{ipa ■ Ua C)q, be an atlas for E, and let faj3 '■ ^p{Ua n Up) (faiUa H C//?) 
denote the corresponding transition functions. Then, Ka/s- Ua r\ Uj3 ^ C* 
given by Ka/sip) = fapi^fsip))'^ is a holomorphic function, that is, 
is an element of the Cech cochain group C^(U, 0*), where U = (Ua) and 
0* denotes the sheaf of non- vanishing holomorphic functions on S. By the 
chain rule, it is actually a cocycle, so it defines an element in H^{U,0*). 
The corresponding holomorphic line bundle K G H^{T,,0*) is called the 
canonical bundle over S. With the notations of Section [21 K is nothing 
but the holomorphic cotangent bundle T*E"'", while K coincides with T*T,~ . 
Hence, the d operator can be seen as a map d: C°°(l) — )• C°^{K), where 1 
denotes the trivial line bundle. 

The set S(S) of spin structures on S can be defined as the set of iso- 
morphism classes of holomorphic line bundles that are square roots of K, 
that is, 

§(S) = {L G 0*) I = K}. 

This is easily seen to be an affine space over H^{T,;7j2). Note that a spin 
structure L is given by a cocycle {Xa/s) G Z^(U,0*) such that A^^ = Hap- 
Then, a spinor ip G C°°{L) can be described by a family of smooth functions 
■00 € C°°(C/«) such that ipaip) = K/^ip) i^fsip) for p G C/« n C/g. Since A«/3 is 
holomorphic, the assignement (ipa) ^ (dipa) defines a map 

called the twisted d operator. An element tp G C°^{L) is a holomorphic 
spinor if it is in the kernel of di- 

Finally, the choice of a hermitian metric on S allows to define an anti- 
linear isomorphism h: C^{L K) — >■ C°°{L). The Dirac operator is the 
self-adjoint operator on C°°(L)©C°°(L) whose restriction to C°°{L) is given 

by 

DL = hodL: C°°{L) C^{L). 
By abuse of language, we shall also call the Dirac operator. 

Let us now give discrete analogs of the objects described above. As ex- 
plained in the previous section, an analog of a Riemann surface is a bipartite 
graph r embedded in a fiat surface S with cone type singularities supported 
at S, inducing a cell decomposition AT of S. Furthermore, in order to define 
d: — )• Q}{W), we assumed that F is isoradially embedded in S, with S 
contained in B U V{T*). 

By definition, Sq := S \ 5 is endowed with an atlas whose transition 
functions are Euclidean isometrics. Therefore, the associated Cech cocycle 
consists of S^-valued constant functions. This defines an element Kq of 
H^{T,Q]S^). Using the long exact sequence for the pair (S,So), one easily 
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checks that Kq is the restriction of a class [k] G H^(Tj; S^) = H^{X; S^) if 
and only if exp{i6x) = 1 for all x E S". We shall therefore assume that all cone 
angles 6x are positive multiples of 27r, and call such a cocycle k G Z^{X; S^) a 
discrete canonical bundle over S. Note that the cohomology class of k is 
uniquely determined by the flat metric on E and the cellular decomposition 
X. Furthermore, such a cocycle k is very easy to compute, as demonstrated 
by the following remarks. 

Remark 3.1. If the flat surface E has trivial holonomy, one can simply choose 
K = 1 as discrete canonical bundle. 

Remark 3.2. It is always possible to represent S as planar polygons P with 
boundary identifications. Furthermore, these polygons can be chosen so 
that r intersects dP transversally, except at possible singularities in S (1 B. 
Define k by 



where ?9 denotes the angle between the sides of dP C C met by the edge 
e. If S is contained in V{T*), this defines completely a natural choice of 
discrete canonical bundle k. If S* fl S is not empty, the partially defined k 
above can be extended to a cocycle yielding a discrete canonical bundle. 

Example 3.3. Let P be the regular 4g gon with boundary identification 
according to the word Yl^j=i0.jbjaj^bj^ . This defines a flat metric on the 
genus g orientable surface T,g with one singularity of angle 2TT{2g — l). Given 
a graph F C meeting dP transversally, the associated canonical bundle 
is given by ^(e) = exp(— ivr^^) for edges of F meeting dP, and K(e) = 1 for 
interior edges. 

Finally, note that the assumption that all cone angles are multiples of 
Itt simply means that E has trivial local holonomy. In such Si CclSG, the 
holonomy defines an element of IIom(7ri(E), S^) = -ff^(E; S^) = H^{X; S^), 
and a representative of this cohomology class is exactly the inverse of a 
discrete canonical bundle. Note also that this assumption rules out the 
2-sphere from our setting. 

Mimicking the continuous case, let us define a discrete spin structure 
on E as any cellular 1-cocycle A G Z^{X;S^) such that = k. (See [5j for 
another notion of discrete spin structure, valid for any cellular decomposi- 
tion of an orientable surface.) Two discrete spin structures will be called 
equivalent if they are cohomologous. The set §{X) of equivalent classes of 
discrete spin structures on E is then given by 



One easily checks that this set admits a freely transitive action of the abelian 
group i7^(E; {±1}). In other words, and using additive notations, §{X) is 




1 if e is contained in the interior of P; 

exp(— if e meets dP transversally. 



S{X) = {[X]eHHx-S')\[Xf = [K]}. 
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an affine H^{T,; Z2)-space. Therefore, there exist (non-canonical) H^{Ti; Z2)- 
equivariant bijections §{X) —?■ S(5]). Furthermore: 

Proposition 3.4. If all cone angles ofT, are odd multiples of 2Tr, then there 
exists a canonical H^{Y!,;Zi2)-equivariant bijection §{X) — )• S(S). 

Proof. Let k £ Z^{X;S^) be a fixed discrete canonical bundle over S. For 
each A G Z^{X; S^) such that = k, we shall now construct a vector field 
Vx on S with zeroes of even index. Such a vector field is well-known to define 
a spin structure, or equivalently - by Johnson's theorem [11] - a quadratic 
form q\ on Hi{T,; Z2). The proof will be completed with the verification that 
two equivalent A's induce identical quadratic forms, and that the assignment 
[A] I— 7- qx is Z2)-equivariant. 

Let A G Z^{X; S^) be given by A(e) = exp{ij3x{e)) with < /3a (e) < 27r, 
where e is an edge of X oriented from the white end to the black end, and 
set (3x{-e) = -(3xie). 

First, replace the cellular decomposition X 
of S by X' , where each singularity b & B n S 
is removed as illustrated opposite. Obviously, 
A induces A' G Z\X'; S^) by setting A'(e) = 1 
for each newly created edge e. Now, fix an arbitrary tangent vector Vx{w) 
at some white vertex w, and extend it to the 1-skeleton F' of X' as follows: 
running along an edge e oriented from the white end to the black end, rotate 
the tangent vector by an angle of 2/3^(6) in the negative direction. (On the 
newly created edges, just extend the vector field without any rotation.) As 
A' is a cocycle and each cone angle is a multiple of 27r, this gives a well- 
defined vector field along F'. Extend it to the whole surface S by the cone 
construction. The resulting vector field Vx has one zero in the center of each 
face of X, and at each b G B S. One easily checks that such a zero is of 
even index if and only if the corresponding cone angle is an odd multiple of 
2tt, which we assumed. 

Following [llj, the quadratic form qx'. -ffi(Il;Z2) Z2 corresponding to 
Vx is determined as follows: for any regular oriented simple closed curve 
C C S \ 5, gA([C']) + 1 is equal to the winding number of the tangential 
vector field along C with respect to the vector field Vx- For an oriented 
simple closed curve C C F, we obtain the following equality modulo 2: 

9A([q) = l + ^(^2/3,(e) + 5;(vr-a.(C))) 

ecC veC 

= l + ^ + ^(E2/3A(e)-E«^(^))' 

eCC v&C 
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c 

where the first sum is over all oriented edges in the ori- ^ av{C) 
ented curve C, and ayiC) is the angle illustrated oppo- ^C^^ 
site. Obviously, equivalent A's induce the same quadratic ^ 
form q\. Finally, given two discrete spin structures Ai,A2, the cohomology 
class of the 1-cocycle A1/A2 G Z^(X;{zbl}) is determined by its value on 
oriented simple closed curves in F. For such a curve C, we have 

(Ai/A2)(C) = exp (iY^{fixM - /3A.(e))) = exp ^m{q^, - (?a.)([C])). 

eCC 

Therefore, the assignement [A] 1-^ q\ is Z2)-equivariant, which con- 

cludes the proof. □ 

Remark 3.5. If the flat surface S has trivial holonomy, then [k] is trivial, so 
the set S(X) is equal to the 2(7-dimensional vector space {Yj^Ia^) . 

Remark 3.6. Let F C S be described via planar polygons as explained in 
Remark [3. 21 and let us assume that the singular set S is contained in F(F*). 
In such a case, a discrete spin structure is given by 



A(e) 




if e is contained in the interior of P\ 
-i'd/2) if e meets dP, 



where exp(— i??/2) denotes one of the square roots of the angle between the 
sides of dP C C met by the edge e. 



Example 3.7. For S as in Example 13.31 above, equivalence classes of spin 
structures correspond to the 2^^ choices of 2g square roots of exp(— ivr^^), 
one for each pair of boundary edges of P. In particular, for the flat torus, 
§>{X) corresponds to the 4 possible choices of 2 square roots of the unity. 



Let us now turn to spinors. Given a spin structure L G S(S), the universal 
covering vr: S — t- S induces the following pullback diagram: 

n 



■K p 




By the lifting property of the covering map H: E L, any spinor S 
C°^{L) induces a section il) £ C°°{E) such that li o if) = o -k, unique up 
to the action of 7ri(E). Since S is contractible (recall that the 2-sphere is 
ruled out by the assumption on the cone angles), the line bundle E is 
trivial. Hence, ijj is really a complex-valued function on S satisfying some 
7ri(S)-periodicity property depending on L. This alternative point of view 
on spinors leads to the following definition. 

Let A € Z^{X] S"^) be a discrete spin structure on S, and let vr: X — > X 
denote the cellular map given by the universal covering of S. Note that 
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Table 2. Discretization dictionary, part 2 



the geometric object 


the discrete analog 


the canonical bundle 
K G H\T.,0*) 


the cohomology class [k] € H^{X; S^ ), 
provided all cone angles {9x}x&s are 
multiples of 2tt 


the afhne 
H\^-Z2)-space S(S) of 
spin structures 


the affine //""^(S; Z2)-space S(X) of 
equivalence classes of square roots of 
K € Z^{X-S^) 


the spinors C°°{L) 
associated to L G S(S) 


the space C(A) C C'^ of discrete spinors 
associated to A € S(^) 




C(\) C 


a hermitian metric on S 


a (normalized) vector field V £ X(W^) 


the twisted d operator 


5a: C(A) — )• r2-^(VK) given by the restriction 
of d to C{\) C C-^ 


the Dirac operator 


Da : ~ C7(A) ^ C7(A) ~ given by ^a 
evaluated along the vector field V 



the bipartite structure on F lifts to a bipartite structure ^(r) = B \JW on 
r = 7r~^(r). The space C(A) of discrete spinors is the set of all -0 G 
such that, for any b,b' & B with 7r(&) = vr(6'), 

V'(6') = AM7M'))V'W, 

where 7f,,b' denotes a path in F from b to 6'. As A is a cocycle and T, is 
simply-connected, this condition does not depend on the choice of such a 
path. Furthermore, equivalent discrete spin structures A ~ A' will yield the 
same space C(A) = C(A'). Note that the choice of any fundamental domain 

P C S for the action of 7ri(S) yields an identification C(A) ~ C . However, 
this identification is not canonical, unless A is trivial. 

Similarly, let us define the space C(A) as the set of all ip G such that 
tp{w') = X{TT{'yw,w')) tpiw) whenever Tr{w) = Tr{w'). Here again, a fundamen- 

tal domain P C T, yields a non-canonical identification C(A) ~ C . 

Finally, and for reasons that will become clear in the next paragraph, the 
role of the hermitian metric on S will be played by a nowhere vanishing 
vector field V G X{W) along the white vertices of F. Furthermore, we shall 
normalize this vector field so that {V^l = 2Area{St{w)) for all w £ W. 

Table [2] summarizes the second part of our dictionary. 
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3.2. The discrete Dirac operators. As above, let F be a bipartite graph 
isoradially embedded in a flat surface S with cone type singularities S C 
B U V{T*), and let us assume that all cone angles are multiples of 27r. Note 
that all these structures lift to the universal cover vr: S — )• S. Indeed, this 
map defines a bipartite graph T isoradially embedded in the flat surface S 
with cone type singularities S C B U V(T*). Let us define the discrete 
twisted d operator associated to A G S(-'^) is the C- linear map 

defined by the restriction of the discrete d operator B: Q^{W) to 

C{X) C C^. 

We need a map r2^(W^) — > discretizing the anti-linear isomorphism 
C°°{L (8) K) — )• C°°{L) induced by a hermitian metric. A discrete hermitian 
metric, that is, a normalized vector field V G induces a very natural 

such map, namely the evaluation at F S X{W), the lift of V to W. Putting 
all the pieces together yields the map D'^: C(A) — > given by 

with the notations of Section 12. 31 One easily checks that the image of D'^ 
is contained in C(A), and that equivalent discrete spin structures A ~ A' 
induce identical maps D'-^ = D'^, : C(A) — >• C(A). Finally, the following 
lemma provides us with a less cumbersome definition of this operator. 

Lemma 3.8. Pick a simply- connected fundamental domain P C S for the 
action of Tri{T,), and let D\: — > be the composition (pp o D'^ o tpp^. 
Then, for a well-chosen representative of [A] G ^(X), 

{Dxi^){w) = ^X{w,b)u{w,b)e''^''^'"'^^i;{b) 
for ip G C'^ and w G W. 

Proof. Fix ip G C'^, w G W, and let w denote the element of n^^{w) in P. 
Then, 

= Y^viw^^e'^'v^^^^^ ^-p\^){b) 

= Yuiwrb)e''^v^^''H{n{r,>;,mAb)^ 

where b' denotes the element of P such that 7r(6') = vr(6). As P is simply- 
connected, there exists a representative A such that A(e) = 1 for any edge e 
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contained in the interior of tt{P). Setting 7r(6) = 6, we get 

what was to be shown. □ 

This discussion motivates the following definition. 

Definition 3.9. Let F be a bipartite graph isoradially embedded in a flat 
surface S with conical singularities S C y(r*) U -B, and all cone angles 
multiples of 2ti. Given any discrete spin structure A, the associated discrete 
Dirac operator is the map D\ : C'^ defined by 

for ^ G and t(; G W. The sum is over all vertices h adjacent to w, v{w,h) 
denote the length of the edge dual to {w,b), and {^yiw^h) is the angle at 
w illustrated in Figure [3l 

A discrete spinor G is discrete holomorphic (with respect to 
A) if Dxip = 0. 

Note that D\ is essentially independant from the choice of the discrete 
hermitian metric V: another choice would yield the matrix QDx, where Q 
is a diagonal matrix with diagonal coefficients in S^. Furthermore, if A and 
A' are equivalent discrete spin structures, then there exist two such matrices 
Q, Q' such that Dx' = QDxQ' ■ 

Remark 3.10. The map Dx'- — ?• defined above is really the discrete 
analog of the restriction of the Dirac operator to C°°(L). The full Dirac 
operator on C°°(L) © C°°{L) being self-adjoint, it would discretize to the 

operator on C^^^^ = © C'^ given by the matrix 

Remark 3.11. We have assumed throughout the paper that no white vertex 
of F is a conical singularity of S. This was crucial in Section [2] in order to 
define the discrete d operator. However, in the present section, we could 
have dropped this condition and defined Dx using any choice of a direction 
at each w € W (for example, given by a perfect matching). All the results 
of the paper, apart from the ones of Section [21 still hold in this slightly more 
general setting. 

3.3. The convergence theorem for spinors. Let us conclude this section 
with the application of the convergence theorem (Theorem 12. 5p to spinors. 

Let F„ be a sequence of graphs embedded in a flat surface S, and let 
A„ G §>{Xn) be discrete spin structures inducing the same spin structure 
L G S(S). (Recall that by Proposition 13.41 there is a canonical equivariant 
bijection §{X) — )■ S(S) provided all cone angles are odd multiples of 27r.) We 
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shall say that a sequence V'n £ C{Xn) C C^" of discrete spinors converges 
to a section ip of the line bundle L ^ S if, for some lift ip : T, C of , 
the following holds: for any sequence x„ £ Bn converging to x S S, i^n{xn) 
converges to ip{x). 

Theorem 3.12. Let Ti be a flat surface with conical singularities supported 
at S whose angles are odd multiples of 27r. Consider a sequence r„ of bi- 
partite graphs isoradially embedded in S with S C V^(r*) U B^, inducing 
cellular decompositions Xn ofTi. Assume that the radii 5n ofV^ converge 
to 0, and that there is some r] > such that all rhombi angles of all these 
Tn 's belong to [rj, n — rj]. Finally, pick a sequence of discrete spin structures 
Xn € Z^{Xn;S^) inducing the same class in H^{Ti;S^), and let L G S(S) 
denote the corresponding spin structure on E. 

Let ipn £ C'(An) be a sequence of discrete spinors converging to a section 
ip of the line bundle L — )• S. If for each n, ipn is discrete holomorphic with 
respect to Xn, then ijj is a holomorphic spinor. 

Proof. By assumption, ipn G C^" are discrete holomorphic functions on S 
converging to ip: S ^ C in the sense of Theorem 12.51 By this result, ip is a 
holomorphic function. Therefore, ip G C°°{L) is a holomorphic spinor. □ 



4. Relation to Kasteleyn matrices and the dimer model 



Recall that a dimer covering, or perfect matching on a finite con- 
nected graph r is a family M of edges of F, called dimers, such that 
each vertex of F is adjacent to exactly one dimer. Any edge weight sys- 
tem v: E{r) [0, oo) induces a probability measure on the set M(r) of 
dimer coverings of F. It is given by 

Z{T,u) 

where i>{M) = HeeA/ '^i^) 

ZiT,u)= KM) 
AfeM(r) 

is the associated partition function. The study of this measure is called 
the dimer model on F. 

The aim of this section is to relate the discrete Dirac operators introduced 
above to some matrices, called Kasteleyn matrices, which provide a standard 
tool for the dimer model on a graph. 

4.1. Kasteleyn flatness. Let F be a finite bipartite graph. Fix a field F 
containing R as a subfield, and let G be a multiplicative subgroup of F* 
containing {±1}. (The examples to keep in mind are G = {±1} C M* and 
G = C C* .) Since each edge of F is endowed with a natural orientation 



26 



David Cimasoni 



(say, from the white vertex to the black one), a map lo: E(T) — >■ G can be 
viewed as a cellular 1-cochain oj G C^{T;G), where co ( o — » • • ) = uj{e) 
and CO ( cv^* — • ) = uj{e)~^. 

Let us order the set B of black vertices of F, as well as the white vertices 
W, and fix a cochain w G C^(r;G). Let iC'^ = iC'^(r,z/) denote the asso- 
ciated weighted bipartite adjacency matrix: This is the (|VF| x |iJ|)-matrix 
with coefficients in ¥ defined by 

e 

the sum being on all edges e of F joining w eW and b E B. 

The goal is now to find cochains oo so that det(i^'^(r, z^)) can be used 
to compute Z(r, i/). Embed F in an oriented closed surface S so that S \ 
F consists of open 2-discs (this is always possible), and let X denote the 
induced cellular decomposition of S. The Kasteleyn curvature of w € 
C^(F; G) at a face / of X is the element of G defined by 

where df denotes the oriented boundary of the oriented face /, and \df\ the 
number of edges in df. This defines a curvature 2-cochain c^^ € C'^{X; G). 
A 1-cochain cu is said to be Kasteleyn fiat (or simply fiat) if c^^ is equal 
to L Finally, we shall say that two cochains oj,uj' G C^(F;G) are gauge 
equivalent (or simply equivalent) if they arc cohomologous, that is, if 
they can be related by iterations of the following transformation: pick a 
vertex of F and multiply all adjacent edge weights by some g & G. Note that 
equivalent cochains uj,uj' have the same curvature, and that the determinant 
of the associated matrices K'^ and K'^ differ by multiplication by an element 
of G. 

Example 4.1. If G is the multiplicative group {±1}, then elements of (F; G) 
are nothing but orientations of the edges of F: an edge e is oriented from the 
white vertex to the black one if and only if a;(e) = -|-1. Furthermore, oj is flat 
if and only if the corresponding orientation satisfies the following condition: 
for each face /, the number of boundary edges oriented from black to white 
has the parity of + 1. This is usually called a Kasteleyn orientation, 
and the associated matrix K'^ is called a Kasteleyn matrix. By abuse 
of language, we shall say that two Kasteleyn matrices are equivalent if the 
corresponding Kasteleyn orientations are. 

Proposition 4.2. There exists aflat G-valued 1-cochain on a bipartite graph 
r C T, if and only if F has an even number of vertices. In this case, the set 
of equivalence classes of such 1-cochains is an H^{T,;G)-torsor, that is: it 
admits a freely transitive action of the abelian group i7^(S;G). 
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Proof. Let V (resp. E, F) denote the number of vertices (resp. edges, faces) 
of X. Given any io € C^{T; G), we have 

n cM) = = (-1)^+^ = (-l)^ 

fcx 

since the Euler characteristic x(^) = V—E+F is even. Therefore, if uj is flat, 
then V is even. Conversely, if V is even, then Il/cx ^i^if) ~ 1- This implies 
that Cu) is a coboundary, that is, there exists a (j) £ C^{X]G) such that 
Cuj = <5</>~^. Consider now the 1-cochain (j)Lo defined by {cl)Lo){e) = cl){e)Lo{e). 
Given any face / of X, we have the following equality in G: 

Since c^^ = 5(j)^^, it follows that c^p^ = 1, that is, (puj is flat. 

Let us now prove the second statement, assuming that there exists a flat 
cochain. Define the action of an element [(p] € G) = H^(X] G) on [w] 

by [0] • [a;] = [</>w]. Since (/> is a cocycle, the equation displayed above shows 
that (f)UJ is flat if and only \i u is. Note also that is gauge equivalent to 
UJ if and only if (/> is a coboundary. Therefore, this action of H^(Ti]G) on 
the set of equivalence classes is well-defined, and free. Finally, given two flat 
systems uj and uj', let (j) denote the 1-cochain defined by (/>(e) = uj'{e)uj{e)~^ . 
Obviously, uj' = (j)UJ, and is a cocycle by the identity displayed above. 
Therefore, the action is freely transitive. □ 

4.2. Computing the dimer partition function. The point of introduc- 
ing flat cochains is that they can be used to compute the partition function 
Z(r, v) of the dimer model, as follows. Note that Z{T, v) is zero unless V 
has the same number of white and black vertices, which we shall assume 
throughout this section. 

Let 23 = {aj} be a set of simple closed curves on S, transverse to F, whose 
classes form a basis of -f^i(S; Z). For each aj G let Gj denote the oriented 
1-cycle in F having aj to its immediate left, and meeting every vertex of F 
adjacent to aj on this side. Let r denote the flat cochain (unique up to 
equivalence, by Proposition 14. 2p such that T{Gj) = (— 1)1*^^1/2+1 £qj. ^ 
Let UJ G C^(F; G) be any flat cochain, and let (p be the unique element in 
H^{Tj; G) such that Lp -[t] = [uj]. Finally, for any e = (ei, . . . , e2g) € Zg^, let 
We denote the flat cochain obtained from uj as follows: multiply uj{e) by — 1 
each time the edge e meets aj with ej = 1. 

Theorem 4.3. For any a G i^i(S;Z), let Z'^{T,v) denote the partial par- 
tition function defined by 

MGM(r) 
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Then, the following equality holds in ¥ up to multiplication by an element 
ofG: 

Taking the flat cochain u = t immediately yields: 
Corollary 4.4. The partition function of the dimer model on T is given by 

^(r,^) = ^ E (-l)^'<^"^^^"""^det(E:^^). 

where = stands for equality in ¥ up to multiplication by an element ofG. □ 

Example 4.5. Assume that the bipartite graph T is planar. In such a case, 
one can take S to be the 2-sphere, so all flat G- valued cochains are equivalent 
by Proposition 14.21 Corollary 14.41 gives the equality 

Z{T,v) = det(i^'^) 

for any flat uj € C^{V;G). The case G = {±1} is the celebrated Kasteleyn 
Theorem |121 [T3t . The mild generalization stated in this example is not 
truely original, as it easily follows from the discussion in Section II of [16j. 

The general formula stated in Theorem 14.31 can seem somewhat cumber- 
some. Therefore, let us illustrate its usefulness before giving the proof. 

Example 4.6. Let F be the quotient field of the group ring Z[Hi{T,;Z)], and 
let G denote the subgroup of F* given by G = ±Hi{T,;Z). If one chooses 
a family of curves 23 = {oi} as above and denotes by € the 
class of ai, then F = Q(ai, . . . ,a2g). Let r € G^(r; {±1}) be as described 
above, and consider the cochain uj £ G^(r; G) given by uj{e) = r(e) Y\ - a"''*^. 
In other words, an edge is multiplied by Oj (resp. a~^) each time it crosses 
ai in the positive (resp. negative) direction. Then, Theorem 14.31 vields the 
following equality in Q(ai, . . . , a2g), up to multiplication by ita™^ ■ ■ ■ a^g": 

Zl{T,^)aT---all^ = 1 j;(-l)^-.-^^^-'-^det(i^-^), 

where Z^{r, u) is the partial partition function given by the contribution of 
ah M G M(r) such that • M = n, for all i. 

Proof of Theorem \4-3\ Consider a itl-valued cochain r, and interpret it as 
an orientation of the edges of T as explained in Example 14.11 One easily 

checks that the equation r(Cj) = (— 1)~2 ^ is equivalent to the following 



Discrete Dirac operators and Kasteleyn matrices 



29 



fact: the number of edges in Cj where r disagrees with a given orientation 
on Cj is odd. By [U Theorem 3.9], we have the following equality in F 

z(r,w) = ^Y1 (-i)^^<^"'^^"""^ det(i^-^(r,w)), 

for any F- valued edge weight system w. Given any cochains a^cf) ^ C^{^] G), 
the equality 

is obvious. Furthermore, if (/> is a cocycle, we shall check shortly that 

where [cj)] G H^{T,;G) = Hom(ffi(S; Z), G) is the cohomology class of (p, 
and Z^(r,z^) is the partial partition function defined in the statement of 
the theorem. Applying the three equalities displayed above to the weight 
system w = with cj) such that (j)T = uj, yields the theorem. 

It remains to check the last equation displayed above. Let {a*} be the 
basis in Z) dual to S = {aj} with respect to the intersection pairing, 

that is, such that Oj • = 5ij. The difference of any two dimer coverings 
M, Mq viewed as elements of Ci{T; Z) is clearly a cycle. Since the expression 
of an arbitrary a G Hi{T,; Z) in the basis {a*} is given by a = Ylii'^i ' 
we get 

This implies the equality 

AfGM(r) MeM(r) i 

aGHi(S;Z) 

which concludes the proof. □ 

4.3. Discrete Dirac operators and Kasteleyn matrices. Now, let us 
turn back to our discrete Dirac operators. As in Section[2l let F be a bipartite 
graph isoradially embedded is a flat surface S with cone type singularities 
5 CBUV(T*). 

Lemma 4.7. Given a nowhere vanishing vector field V along W , let LOy G 
C'^iT] S"^) be the cochain defined by ujv{^) = exp(i'!?y (u), 6)) as illustrated in 
Figure\^ Then, the equivalence class of coy does not depend on V. Further- 
more, its Kasteleyn curvature is given by 

Cojvif) = -exp(i6l//2). 
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where Of denotes the angle of the conical singularity in the face f . 

Proof. The first statement is obvious, so let us fix a nowhere vanishing vector 
field V G X(VF) and consider the associated cochain wy. Given a face / 
of r C S, let wi,hi,W2,h2, ■ ■ ■ ,Wm,bm denote the vertices in df cyclically 
ordered. Then, the Kasteleyn curvature of ajy the face / is given by 

= {-l)'"+^u)v{wi,bi)u:v{bi,W2)ujv{w2,b2) ■ ■■ujv{bm,wi) 

- ( I)'" ^V{wi,bi)ujv{w2, ^2) • • • tdv{Wm, bm) 

UJviwi,bm)uJviw2, bl)--- UJviWm, &m-l) 
m 

= -exp (i^(7r - Q^^(a/))), 
i=i 

where aw.{df) = '&v{wj,bj-i) — '&v{'Wj,bj) and 60 = bm- This angle a^idf) 
is simply the angle made by the oriented curve df at the vertex w, as 
illustrated in Figure [H An easy application of Gauss-Bonnet shows that the 
angle 6*/ of the conical singularity x/ in / is equal to J2vedf('^ ~ '^vidf))- 
Hence, it remains to check that 

^ a,{df)- «»(5/) = 0- 

beBndf wewndf 

This is where the isoradiality comes into play. By definition, there is a 
local isometry from the pointed face f \{xf} to the pointed plane C* such 
that all vertices in df are mapped to a circle in C*. Now, observe that the 
alternating sum of angles displayed above does not change if one moves a 
vertex along the circle keeping all other vertices fixed. Since the equality 
above holds when all angles are equal (to vr — Of /2m), this concludes the 
proof. □ 




V 



Figure 4. The angle made by the oriented curve C at the 
vertex v. 

Let us now state the main result of this section. 

Theorem 4.8. Let T, be a compact oriented fiat surface of genus g with 
conical singularities supported at S and cone angles multiples of 2tt. Fix a 
graph T with bipartite structure V{T) = B Li W , isoradially embedded in S 
so that S C BU V{T*). For an edge e ofT, let v{e) denote the length of the 
dual edge. Finally, let D\ : — >■ denote the discrete Dirac operator 
associated to the discrete spin structure X. 
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There exist "2}^ non- equivalent discrete spin structures such that the corre- 
sponding discrete Dirac operators {Dx}x give non- equivalent Kasteleyn 
matrices of the weighted graph (r, i/), if and only if the following conditions 
hold: 

(i) each conical singularity in V{T*) has angle an odd multiple of2iT; 
(a) for some (or equivalently, for any) choice of oriented simple closed 
curves {Cj} in T representing a basis of Hi{T,;'Z), 

b€BnCj wGWnCj 

is a multiple of 27: for all j , where a^,(C) denotes the angle made by 
the oriented curve C at the vertex v as illustrated in Figure [7} 

Proof. Fix a discrete spin structure A € Z^{X] S"^), a normalized vector 
field V E ?i{W), and let D\: — > be the corresponding discrete 
Dirac operator. By definition, the coefficient of Dx corresponding to vertices 
w G W and 6 € -B is equal to 



Dx{w,b) 



i>{e) ujv{e)X{e) if w and b are joined by an edge e; 
ii w and b are not adjacent, 

with uJv ^ in Lemma 14.71 In other words, Dx is the adjacency matrix of 
the weighted bipartite graph (r,z^), twisted by the cochain ojx ■= wyA € 
C^{T; S^). The goal is now to check that ux is gauge equivalent to a Kaste- 
leyn orientation (that is, to a ±1- valued flat cochain) if and only if Condi- 
tions (i) and (ii) hold. This clearly implies the theorem, as non-equivalent 
discrete spin structures yield non-equivalent Kasteleyn orientations. 

By Lemma [4. 71 LOy is flat if and only if Condition (i) is satisfied. Since A is 
a cocycle, 6X = 1 and the same statement holds for ojx- By Proposition 14.21 
the set of equivalence classes of such S^-valued flat cochains is an H^{T,; S^)- 
torsor. Therefore, ujx is equivalent to a Kasteleyn orientation if and only 
if, for any Kasteleyn orientation ojq, the cocycle (f) := ojq^ujx represents a 
class in H^{T.;{±1}) = Hom(Fi(S; Z), {±1}). This holds if and only if 
0(C) € {il} for any 1-cycle C in F, or equivalently, for any 1-cycle in F 
part of a collection representing a basis of /7i(S;Z). Since ujq = 1, this 
translates into the equalities 

1 = cj){Cf = ujxiCf = ujv{C)'^X{C)'^ = LOviCfuiC). 

As mentioned in Section [3.H k{C) is the inverse of the holonomy along the 
1-cycle C. Therefore, 

k{C) = hol{C)-^ = exp (^-i (vr-a„(C))) = exp (^i ^ a^(C))^ 

v&v{r)nc v£V{r)nc 

since C is of even length. Furthermore, the definition of LOy implies that 
(jJviC) = exp(— i ai„(C)), as in the proof of Lemma [4.71 This 
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yields the equation 

1 = exp (^i ^ ab{C) - i ^ au,(C)^ , 
beBnc wewnc 

obviously equivalent to Condition (ii). □ 



Consider F C S as above, and satisfying both conditions of Theorem [ 
Let 23 = {cxj} be a set of simple closed curves on S, transverse to T, whose 
classes form a basis of Hi{T,; Z). For each aj G 23, let Cj denote the oriented 
1-cycle in F having aj to its immediate left, and meeting every vertex of F 
adjacent to aj on this side. By the conditions of Theorem 14. 8| any discrete 
spin structure A satisfies the equations 



W A(C,)=exp(i Yl 



^+1 



up to a sign. Let us pick the discrete spin structure Aq such that the equality 
above holds for all j. For any e = (ei,...,e2g) € Zg^, let Ag denote the 
discrete spin structure obtained from Aq as follows: 

A,(e) = (-l)^.^^(^-"^)Ao(e). 



Theorem 4.9. If T C S satisfies the conditions of Theorem \4-^ then the 
partition function for the dimer model on (F, v) is given by 

^(r-^) = E (-l)^-^"'^^°'-"^ det(DAj 

Proof. By Condition (i), the /S^-valued cochain Xqujv is Kasteleyn flat. By 
Condition (m), it is gauge equivalent to a { it 1}- valued cocycle r. Finally, 
Equation (*) is equivalent to T{Cj) = (— The theorem now follows 
from Corollary Sai for G = S"! C F* = C*. □ 

Remark 4.10. More generally, let us assume that F C S only satisfies the 
first condition of Theorem 14. 8( and let A be any discrete spin structure. 
Then, Theorem 14.31 gives the equality 

where ip G H'^iT,;!?) is such that 99 • [r] = [Atjy]. 

Spin structures on a closed orientable surface S can be identified with 
quadratic forms, that is, with Z2-valued maps on iJi(E;Z2) such that 
q{x + y) = q{x) + q{y) + x ■ y for all in f/'i(S;Z2). More precisely, 
Johnson [11] gave an explicit Z2)-equivariant bijection S(S) 2(5^) 

between the corresponding affine Z2)-spaces. The Arf invariant of 
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a spin structure is then defined as the Arf invariant of the corresponding 
quadratic form q, that is, the mod 2 integer Arf{q) G Z2 such that 



V|i^i(S;Z2)| 



If all cone angles of S are odd multiples of 2tt, then there exists a canonical 
equivariant bijection §{X) — S(S) (recall Proposition I3.4p . In such a case, 
it makes sense to talk about the Arf invariant Arf (A) of a discrete spin 
structure A € S(-^)- 

As above, let {aj} be a set of simple closed curves on S, transverse to 
r, defining a basis of -ffi(S;Z), and let Cj denote the oriented cycle in T 
having aj to its immediate left. By Condition (ii), the number 

9o(ai) = ^( a^{Cj)- ab{Cj 

wewnCj b€BnCj 

is an integer. Furthermore, one easily checks that its parity changes each 
time aj moves across one vertex. Therefore, the Oj's can be chosen so that 
all (7o(oj)'s are even. 

This leads to the following version of the Pfaffian formula, assuming the 
notations preceding Theorem 14.91 



Theorem 4.11. Let T C T, be as in the statement of Theorem \4.8\ with all 
cone angles ofT, odd multiples o/27r. Then, the partition function for the 
dimer model on (F, u) is given by 
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Proof. We saw in the proof of Proposition 13.41 that the quadratic form 
qx- Z2) —7-^2 corresponding to a class [A] € S(-'^) via the equivariant 

bijection S{X) — > S(S) Q(5^) is determined by the following condition: 
for any oriented simple closed curve C C F, 

qxiiC]) = 1 + ^ + ^ ( E 2/3A(e) - J2 «-(^)) ' 

eCC vac 

where < f^\{e) < 2n is such that A(e) = exp(i/3A(e)). In particular, 
if F C S satisfies Condition (ii) of Theorem 14.81 then the discrete spin 
structure Aq defined by Equation (*) corresponds to the quadratic form qo 
determined by the equalities 

qo{a,) = qoiiCj]) = ^[ ^ a^Cj) - Yl ""(^j))' 

w&wnCj beBnCj 
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By construction, is obtained from Aq by action of the Poincare dual to 
Ae = X^j- ejUj S Hi{T,;Z2). Therefore, by [S, Lemma 1], 

Arf(Ae) + Arf(Ao) = Q'o(Ae) = '^ejqoiaj) + '^eiejai -aj. 

As we have chosen the a^'s so that ah qo^ajYs are even, the formula now 
follows from Theorem 14.91 □ 

Remark 4.12. As mentioned in Remark l3.11l it is not necessary in the present 
section to assume that the sets S and W are disjoint. All the results of this 
section still hold in the slightly more general setting where S is contained 

in y(r) uy(r*). 

4.4. Examples. We conclude this article with a discussion of several special 
cases, and the following result: the Dirac operators on any closed Riemann 
surface of positive genus can be approximated by Kasteleyn matrices. 

The planar case. Let P be a planar isoradial bipartite graph whose associ- 
ated rhombic lattice tiles a simply-connected domain S of the plane. In this 
case, the unique spin structure on S being trivial, the associated discrete 
Dirac operator D : is simply given by 

where the angle '&v{w-, b) can be measured with respect to a constant vector 
field V. With the notations of Figure O this yields the equality 

m 

{D-ip){w) = i ^{xj-i - Xj)f{bj), 
i=i 

which is exactly the discrete Dirac operator introduced by Kenyon |15j in 
this special case. The conditions of Theorem 14.81 being trivially satisfied, D 
is (conjugate to) a Kasteleyn matrix for the dimer model on (P, z^), and the 
associated partition function is given by 

Z{r,u) = |det(D)|. 

Thus, in the planar case, we recover Theorem 10.1 of |15] . 

The genus one case. Let P be a planar isoradial bipartite graph, invariant 
under the action of the lattice A = Zl © Zr C C for some r G H. Fix a 
quadrilateral fundamental domain P C C for this action with P intersecting 
dP transversally, and let P C S = C/A be the corresponding toric graph. 
One of the spin structures on S being trivial, the associated discrete Dirac 
operator D : — ?> is again given by 

{Di>){w) = "^y{w,b)e'^^^'^^^H{b), 
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where V can be chosen to be a constant vector field. The three other discrete 
Dirac operators are obtained from D by multiplying the corresponding coef- 
ficient by —1 whenever an edge crosses the horizontal boundary components 
of P (this gives -Di^o)) the vertical ones (-Do,i), or any boundary component 
(Di_i). These 4 matrices are Kasteleyn matrices if and only if F satisfies 
Condition (ii) in Theorem 14.81 in which case Z(T, u) can be written as an 
alternating sum of the determinant of these matrices. 

Example 4.13. Consider the graph illustrated below. We have inserted next 
to each edge the corresponding coefficient of the matrix D. (The graph is 
normalized so that the sides of P have length 3, V is chosen to be the vertical 
upward direction, and oj stands for exp(27ri/3)). 





1 




|l 




1 / 
















In this example, Condition {ii) in Theorem 14. 8 1 is satisfied. Furthermore, 
one easily checks that the trivial discrete spin structure satisfies Equation 
(*), where ai and Q2 are chosen to be the sides of P. Therefore, Theorem 14. 9 
gives the equality 



#M(F) 



Z{T 

1 
2 



1) = i |det(Z)) + det(Z)i,o) 



+ 28 + 28 - (-28)1 = 42. 



+ det(i:'o,i) - det(i:» 



An example of genus 2. Consider the flat surface S of genus 2 given by 
an octagon, where all pairs of opposite sides are identified. Embed a square 
lattice F in S as illustrated below. 

The flat surface S has a single singularity, which lies in V{T*), and has 
angle Gvr. Therefore, this example satisfies the first condition of Theorem l4.81 
One easily checks that it also satifies the second condition, so the 16 discrete 
Dirac operators are Kasteleyn matrices. Since S has trivial holonomy, one 
discrete spin structure can be chosen to be trivial. Furthermore, one can fix 
a constant direction y on S (say, to the right). The corresponding discrete 
Dirac operator D : — t- is simply given by 

{DtP){w) = i^Pih) - ^Pib3)) + i(V(&2) - ^^(^4)), 
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where bi is the black vertex to the right of w, 62 above, 63 to the left, and 64 
below. Using the procedure described before Theorem l4.9l it is now a trivial 
matter to write the number of dimer coverings of T as some alternating sum 
of determinants of these 16 discrete Dirac operators. 

The example above only discretizes the Dirac operators on one specific 
Riemann surface of genus 2. Can one find examples for any closed Riemann 
surface? Obviously not for the Riemann sphere, as all cone angles are as- 
sumed to be positive multiples of 27r. However, this turns out to be the only 
exception, as demonstrated by the following theorem. 

Theorem 4.14. For any closed Riemann surface of positive genus, there ex- 
ist a flat surface S with cone type singularities inducing this complex struc- 
ture, and an isoradially embedded bipartite graph F C S, with arbitrarily 
small radius, satisfying all the hypothesis and conditions of Theorem \4.8\ 

Proof. The building block of our construction will be the rhombus consisting 
of two equilateral triangles glued along one of their sides. Given positive 
integers n and m, let R{m, n) denote m rows of n such rhombi stacked in 
the following way. (This picture represents i?(2,8).) 




Let r be the associated bipartite isoradial planar graph, where the bottom 
left corner of R{m, n) is a black vertex of F. If n is even, the identification 
of the two vertical sides of R{m, n) will preserve the bipartite structure of 
r. If m is even, one can also identify the horizontal sides, possibly with a 
shift. This allows to realize any torus C/Z + Zr with r in some dense subset 
of H. To obtain all tori, continuously deform one or two rows of rhombi as 
illustrated below. 

The deformation of two rows changes the imaginary part of r, while the 
deformation of a single row changes both the imaginary and the real parts. 
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Therefore, a suitable combination of these transformations allows to con- 
struct all tori. These examples are flat tori with no singularity, so they 
trivially satisfy the first condition in Theorem 14. 8[ Furthermore, one easily 
checks that Condition {ii) is also satisfied, provided n and m are divisible 
by 6. (Note that the deformations above do not affect these conditions.) 

Let us now consider a fixed positive even integer n. Given three positive 
integers mi,m2,m3, glue the corresponding rectangles R{mi,n), R{m2,n) 
and R{m3,n) along their bottom side to an equilateral triangle, itself tiled 
by rhombi, as illustrated below. 



Identifying the opposite remaining sides of each rectangle yields a flat pair 
of pants with a single singularity of angle 47r. By varying the m^-'s and using 
the deformation along two rows described above, one can realize any complex 
structure on the pair of pants. Finally, gluing 2g — 2 such pairs of pants along 
their boundaries (with a possible shift and a possible deformation yielding 
a twist), one can realize any Riemann surface of genus g >2. 

To each rhombus, associate the portion of a bipartite graph F given by 

•^■•■..^ J • . (This is just to avoid cumbersome considerations about gluing 

bipartite structures.) The singularities of angle An are located at black 
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vertices of T, so Condition (i) is satisfied. One easily checks that Condition 
(a) is always satisfied for cycles coming from boundary components of the 
pairs of pants. Finally, by chosing wisely the parity of the m^'s, one can 
ensure that Condition (ii) also holds for the cycles passing through several 
pairs of pants. □ 
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